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Abstract 

We study the radial expansion of cylindrical tubes in a hot QGP. These tubes are treated as 
perturbations in the energy density of the system which is formed in heavy ion collisions at RHIC 
and LHC. We start from the equations of relativistic hydrodynamics in two spatial dimensions and 
cylindrical symmetry and perform an expansion of these equations in a small parameter, conserving 
the nonlinearity of the hydro dynamical formalism. We consider both ideal and viscous fluids and 
the latter are studied with a relativistic Navier-Stokes equation. We use the equation of state 
of the MIT bag model. In the case of ideal fluids we obtain a breaking wave equation for the 
energy density fluctuation, which is then solved numerically. We also show that, under certain 
assumptions, perturbations in a relativistic viscous fluid are governed by the Burgers equation. We 
estimate the typical expansion time of the tubes. 
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I. INTRODUCTION 



The study of the initial stage of relativistic heavy ion collisions has experienced a fast 
progress in recent years. One of the most interesting findings in this study, supported 
both by theoretical works and by the analysis of experimental data, is that in the early 
times of these collisions color flux tubes are formed. Although color flux tubes are familiar 
objects in hadron physics it is not obvious that they should be formed in high energy heavy 
ion collisions, where projectile and target can be regarded as bunches of partons without 
any strong clustering neither in configuration nor in color space. Flux tubes, sometimes 
called strings, appear in lattice QCD calculations as field configurations between static 
heavy charges. They appear also in phenomenological models of high energy soft hadronic 
scattering. In the Lund model, for example, when two high energy protons collide with low 
momentum transfer, they cross each other and, due to gluon exchange, a color rearrangement 
takes place with the subsequent formation of strings, which stretch and decay, producing 
particles. 

At high energies, colliding nuclei are dense systems in which the standard linear evolu- 
tion equations (such as the DGLAP equations) must be replaced by others, which include 
nonlinear effects in the evolution. The theoretical description of these dense systems evolved 
into the theory of the Color Glass Condensate (CGC). In this formalism the dense gluonic 
matter is treated in a semi-classical approximation. In this approach, it has been shown in 
l|-|4| that there are solutions of the classical Yang-Mills equations in which the lines of the 
chromo-electric and chromo-magnetic fields are all parallel to the collision axis and these 
fields form color flux tubes in the longitudinal (z) direction. 

The interpretation of RHIC and LHC data also suggests that the system reaches thermal 
equilibrium, forming a thermalized quark-gluon plasma (QGP), very soon after the collision. 
At first sight this would imply that the flux tubes disappear and the quark-gluon matter 
becomes reasonably homogeneous when the hydrodynamical expansion starts. However 
detailed hydrodynamical studies j^, f| strongly suggest that some experimental features 
observed at RHIC and LHC can be understood if we assume that these tubes survive the 
thermalization stage and form "tubular" structures that persist for some time during the 
hydrodynamical expansion. More specifically, the data show the existence of structures in 
the two-particle correlations plotted as function of the pseudorapidity difference At] and the 
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angular spacing A<p. In (5, 6] it has been argued that these structures may have a common 
hydro dynamic origin: the combined effect of longitudinal high energy density tubes (leftover 
from initial particle collisions) and transverse expansion. 

The tubular structures described above, which are nearly uniform in the longitudinal 
direction, may be considered as cylindrical perturbations in the energy density upon a con- 
tinuous background as depicted in Fig. 1. The propagation of perturbations on the top of 



a QGP background has been investigated in several works 



7H10l|. In most of these works 



[3, |8[ a linearized version of hydrodynamics is employed. We have tried to keep the nonlin- 
ear terms in the equations which describe the evolution of the perturbations {9, loj]. This 
extends the validity of our formalism to perturbations which are not so small. 

/ 




FIG. 1: Tubular perturbations on a QGP background. The inner cylinder of radius I represents 
a tube of energy density higher than the background, shown as a cylindrical fireball of radius L. 
The perturbation expands radially. 

In this work we try to answer the question: how fast do the tubes expand in the QGP ? 
In order to obtain the answer we write the hydrodynamical equations for the propagation of 
cylindrical perturbations along the radial direction (see Fig. 1), solve them numerically, and 
estimate what is the time needed for a tube of initial radius of the order of 1 fm to grow and 
reach the typical radius of the system formed in heavy ion collisions, which is of the order 
of 7 fm. If the tube expansion time were much shorter than the lifetime of the fireball, then 
the tube would be very rapidly incorporated in the fireball and it would produce no visible 
effect in the final state particle correlation measurements. 
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We also investigate the effects of viscosity on the expansion of the tubes. It is well 
known that the relativistic version of the Navier-Stokes equation does not constitute a 

n 

causal theory. We strongly recommend the reading of [11] to understand the subject in 
details. Besides these conceptual issues of stability and causality we perform a Navier-Stokes 
approach without worrying about microscopic time scales, due to the nonlinear expansion 
as seen in the appendix. A future and more complete version of this work in relativistic 
viscous hydrodynamics is in progress with the use of Miiller-Israel-Stewart theory, which is 
a causal theory [llf . 

Due to dissipation, viscosity damps the perturbations, which are then more easily mixed 
with the background fluid, loosing their influence on final state particle correlations. 

In contrast to other studies of perturbations in fluids, we do not neglect the nonlinear 
terms in the hydrodynamical equations. 

This text is organized as follows. In the next section we review the basic formulas of 
relativistic hydrodynamics. In section III we review the equation of state (EOS) of the MIT 
bag model. In section IV we derive the equation which describes the evolution of the tube. 
In section V we solve this equation numerically and present some conclusions. 



II. RELATIVISTIC FLUID DYNAMICS 

Pedagogical texts on relativistic hydrodynamics can be found in [ijllisj]. Approximation 



schemes which conserve nonlinearities can be found in [14j and their application to the study 
of nonlinear waves in cold and warm nuclear matter can be found in [15l-ll8| and references 
therein. In this section we briefly review the basic equations (throughout this work we use 
c = 1, h — 1 and the Boltzmann constant is taken to be one, i.e., = 1). 

For simplicity we start our discussion considering two coaxial cylinders. The inner and 
narrower cylinder represents the flux tube which is a perturbation in energy density e. The 
outer and larger cylinder represents the fireball with a uniform energy density 6q (eq < e). 
We will study the expansion of the flux tube in the center of mass system of the fireball. It 
is then natural to chose spatial cylindrical coordinates (z,r,<f)). 

The velocity four-vector u u is defined as u° = 7, u = 7??, where 7 is the Lorentz factor 
7 = (1 — t> 2 ) -1 / 2 and thus u u u u = 1. The velocity field of matter is given by v = v(t, r, z, <p). 
Because of the azimuthal symmetry we do not have components along the direction and 
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consequently no terms involving d/d<fi will survive in what follows. 



A. Ideal fluid 

The energy momentum tensor is given by: 

T,„, = {e + p)u M u„ - pg^ (1) 

where e is the energy density and p the pressure. Energy-momentum conservation is given 
by: 

0„V = (2) 

The projection of (J2]) on a direction perpendicular to yields the relativistic version of 
Euler equation: 

The relativistic version of the continuity equation for the entropy density comes from the 
projection of (^j on the direction of u v : 

(e+ ?) )a/ + /9 )1 e = (4) 

We next recall the Gibbs relation: 

e + p = p B p B + Ts (5) 

and the first law of thermodynamics: 

de = Tds + fisdpB (6) 

In the central rapidity region of heavy ion collisions we expect to find hot QGP with zero 
net baryon number and hence ps = and T ^ 0. Using dps = in ([6]) and inserting ([6]) 
and (j3j) into (jl]) we find: 

Ts(dX) + Tu^s) = 

and finally: 



d v {su v ) = 

as expected for a perfect fluid. This expression can be rewritten as: 



(7) 



B. Viscous fluid 



In order to take the viscosity into account, we add the viscous stress tensor IF" to the 
ideal fluid energy-momentum tensor: 



HQ- 



With this new definition 



where T( ) is the ideal fluid energy-momentum tensor 
of the energy-momentum tensor, Eq. ([2]) remains valid. We will consider 11] a system 
without conserved charges (or at zero chemical potential). As in the case of the ideal fluid, 
we take the appropriate projections of ([2]), which are parallel (u u d^T^ u ) and perpendicular 
(A^d^T^) to the fluid velocity obtaining: 



u v d^ v = De + (e + p)d^ + u v d^W v = 



(10) 



and 



A^T^ = (e + p)Du a - V> + A^IF" = 



where D = u^dn and A^ v = — u^u u . The viscous tensor is given by 



IF" = r]VK v) + (A^V a u a 



where U] 



with 



and 



A {fl B v) = -(A^B V + A V B^) 



V Q = A^du 



(11) 



(12) 



(13) 



(14) 



(15) 



Combining Eqs. f fTOj) and (II ip we can obtain the relativistic version of the Navier-Stokes 



equation. In a compact form it may be found in 
to write it in the long form: 

d 



ljj . For our purposes it is more convenient 



e + p)Y ( + v • V )v + v^- + Vp 



dt 



dt 



d 



d'y - 
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i + v.( T ^) 



\dt 



C- 3^)V 



~dt 



+ V ■ (t#) 



(9t -» 



(16) 



With the help of Eqs. (fTOl) and (llip and using thermodynamical relations we obtain 111 ] : 



For our purposes we shall rewrite it as: 

ds -* -* -* _ r] ( dj 

7^ + 7Vs • v + s — + sV7 • v + 7sV • v = -— I 



(17) 



^ 9 r ^ 
V7- — (7^) 



c?7 - ^ - J 2 
— + 7V • v + V7 ■ V 
at 



(18) 



which is the relativistic version of the continuity equation for the entropy density s. In the 
case of an ideal fluid (rj = £ = 0) we recover the entropy density conservation: 

ds ^^c?7 - _=> _ 

7— + 7 v s • f + s— — h SV7 ■ v + 7sV • v — (J 
at at 

III. EQUATION OF STATE 



From the thermodynamics of the MIT bag model we have 9] : 



87r 2 fi r°° 

3(p + B)=e-B= — T* + - 2 j o dkk*[nj: + n- h 



and 



with the speed of sound c s , given by: 



1 4 « 
p = -e o 

y 3 3 



<9p 1 



(19) 



(20) 



(21) 



de 3 

Since p# = 0, the chemical potential is zero(/i = 0) and the distribution functions are the 
same for quarks and anti-quarks: = ng = 1/(1 + e fc//T ). Therefore: 

37 



3(p + B) = e - B = ^ttT* 



(22) 
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Solving the first identity for the pressure and using the relation s = (dp/dT)v we arrive at: 

s = J* ( - B + -n 2 T 4 ) = 4-tt 2 T 3 (23) 
dT\ 90 / 90 v ; 

The bag constant is related to the critical temperature, T c , of the quark-hadron phase 
transition. During the phase transition the pressure remains constant and (122]) reduces to: 



B 



37 Q T, 



7T 2 — const 



30 3 



37 
30' 



-7T 



30 
377 2 



-const 



and we can define 



and consequently: 



Tn 4 



30 



3 37tt 2 



const 



(24) 



The bag constant, 5, is chosen to be i3 1//4 = 170MeV and this corresponds to Tg = 91 MeV. 
Inserting (I24j) into the second identity of (J22J) we find the following expression for e(T): 

37 



(25) 



Solving the second identity in (I22p for the temperature, we find: 



T 



30 
37^ 



(s-B) 



1/4 



(26) 



which substituted in ( 1231) yields: 



30 
37^2 



13/4 



(e-B) 



From fl22l) we have: 



and from (1211: 



Vp = -Ve 



148 2 . 



and 



dp 1 9e 
dt~3dt 



(27) 



(28) 



(29) 



IV. THE WAVE EQUATION 



In this section we combine the results of the two previous sections and derive the dif- 
ferential equations which govern the evolution of cylindrical perturbations in the energy 



8 



density. We start writing the energy density and the components of the fluid velocity in a 
dimensionless form: 

(30) 



and 



V 

v = — 
c x 



Vz 



(31) 
(32) 



A. Ideal fluid 

After the use of relations (|2"81) to (13"2"j) the components of the Euler equation (J3]) along the 
r and z directions become: 



and 



c 


dv r 
'^t 


2 A dv r 

+ C s V r — 4 

or 


c< 


dv z 
dt 


2 - ^2 

+ c s v r — 4 
or 


d22]) 


, m 


'} and (USD to 


- c s 


2 v 2 ) 


f / 45£ \ " 


|^74vr 2 T 4 / _ 



2 - 

dz 



2 * 



15(c s 2 t) 2 — 1)^0 



747r 2 T 4 V^ + C ^ r at 



(92; 



15(c s 2 ?) 2 — l)e /<9£ * de 

WT 1 lai + Cs ^¥ 



de ^ de A 9s 
— + c s w r — + c s f z — 
at dr dz 



dv\ 



c s v r dv r 

H h c s — + c s 

r dr dz 



,, A <9tV o - ^2 S - 2 <^V 

dt dt dr 



3 ~ 2 



(33) 



(34) 



(35) 



(M 2 3 A A dv r 

~~ ^ 4~ C s V z V r —r i . .s •■- , 

dr dz dz 

Now we combine (1331 . ( )34l) and (|3"5T) to find the wave equation. To this end we perform a 



change of variables in ( 133]) , ( 134]) and ([35 
by the reductive perturbation method 
coordinates [191 ] : 



agoing from the (r, 2, t) space to the (i?, Z, T) space 



19H2l| , through the introduction of the "stretched" 

L - ' " L " ' " L - (36) 

where L is a characteristic length scale of the problem (typically the radius of a heavy ion) 
and a is a small expansion parameter. We next perform the following expansions 19M21|: 



a 1 ' 2 

R = (r - Cgt) , 



Z=°z 
L 



T = -^Cst 



e = 1 + os\ + er 2 £ 2 + cr 3 e 3 + ■ ■■ 



(37) 
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v r = av ri + o 2 v r2 + a 3 v T3 + . . . (38) 

v z = a 3/2 v Zl + a 5/2 v Z2 + a 7/2 v Z3 + ... (39) 

After the use of f[3"6"j) . (I3T1) . f[3"51) and f[3"9"j) . the Euler and continuity equations can be written 
as series in powers of a. We will consider terms only up to the order a 2 . It is then possible 
to reorganize the series in powers of a, a 3 / 2 and a 2 . After a little algebra we find: 

45e \ dei dv ri 1 



(7< 



74tt 2 T 4 / dR OR 



2 J / 45e \ de 2 dv r2 dv ri dv ri ( lhe Q \ de x \ 

M( 45£ ° )^-^Lo (4i) 

1l 747r 2 T 4 i dZ dR K ' 



and 



a< - 



( 45e \de x 
yi^T^jdR ' dR 



dv ri 



2 [ _ ( 45£ \ dE2 ( 45£ \ C>£i / 456Q \ gfi 
\ 1 v 74tt 2 T 4 J <9i? ^74tt 2 T 4 / <9T ^ 74vr 2 T 4 ) Vri OR 

v T1 dv r2 dv Zl _ v ri dv ri \ _ 

T dR dZ 3 dR J ~ U 1 j 

In the above equations each bracket must vanish independently and therefore we obtain a 

set of relations. From the terms of order a in the last two equations we find: 

/ 456 \ dei = dv^ 

\WT^JdR dR 1 } 

which, after the integration over R and taking the integration constant equal to zero, yields: 

^ = ysk) £i (44) 

From the terms of order a 3 ! 2 we have: 



/ 45e \de x 



dv Zl 



V74tt 2 T 4 J dZ dR 
which, after the derivation with respect to Z, becomes: 

( 45 £n \d 2 e x d 2 v 



\n^T±)dZ* dZdR 
From the terms of order a 2 we obtain: 



(45) 



(46) 



dv 



H _ ( 45g ° \ 9e 2 _ du^ ^ dv^ _ ( _L5£o_\ de x 



+ v n^~ WT^i W^b ( 47 ) 



dR \747r 2 T 4 JdR dT ri dR \ 74vr 2 T 4 J ri dR 
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and 



dv r2 
~dR 



( 45e \de 2 ( 45e \de x ( 45e \ dei 
\ 74vr 2 T 4 ) dR ~ \ 74tt 2 T 4 ) dT \ 74tt 2 T 4 ) Vri dR 



i ''Vi dv ri 



T dZ ' 3 dR ^ 
Identifying (j47J) with ( )48l) . using ( 1441) . derivating the resulting equation with respect to i? 
and using (|46|) . we obtain 



5_|5£i + [gf 45 gQ ^ 



3 I 74tt 2 T 4 I 



dei Si 
£l 2T 



+ 



2<9Z 2 







(49) 



Returning now to the (r, t) space we find: 



d_ (dei 



_ -i c dii 
dr \ dt s dr 



° s 51 dr + 2t + 2 dz 2 







(50) 



2/ 45£ V 

3 V 747r 2 r 4 y 

where £ x = cr^i is a small perturbation on the background energy density e . We can 
rewrite the above equation with the coefficients depending only on temperatures and on the 
sound velocity. Using (|2"51) and calling e = e(T = T ), where T is the temperature of the 
background, we then find: 

(51) 



2/ 45e c 



1+ I 



3 V74tt 2 T 4 

where To > Tg. Substituting ( IBTj) into ( 1501) we find the final form of the wave equation: 

}2 



_^_/^£i \- c — | ° s 
<9r [ <9t s (9r 2 



A Sei £i 1 c s g 2 £i 
£l 9r 2t 2 <9z 2 







(52) 



B. Viscous fluid 



Using again the relations ( 128|) to (I32|) in ( 1161) . performing the same change of variables 
( |36|) . performing the same expansions f )37|) . ( 138]) and (139]) . organizing the several terms in 
powers of cr and obtaining the corresponding identities and returning to the (r, z, t) we arrive 
at the analogue of fl50|) for a viscous fluid: 

2/ 45£ 



d (dii dii 
— { — + c s — - + 
dr [ dt s dr 



3 \, 747^^ 



c s e r 



dii ii ( 



45 



dr 2t \ 74tt 2 T 



' 4 \ d 2 iA c s d 2 ii 

(53) 

As expected, the above equation reduces to the corresponding equation for ideal fluids f )50|) 
in the limit 77 = ( = 0. Since the derivation of the equation is very similar to the sequence 
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of steps that led to (1501) . we omitted all the details. However the interested reader can find 
some more details in the Appendix. 

With the help of ( 12 3 p we can rewrite the viscosity term as a function of the dimensionless 



ratios rj/s and £/s, which are well studied in the literature 
( 15"3"|) becomes, after using (lo"Tj) and 



22 



23]. So the wave equation 



1+ I 



£l * 2i T L 3 s J <9r 2 



c s <9 2 £i 

2 a? 2 



(54) 



C. Effect of the background expansion 



So far we have considered the motion of a perturbation on a static background. In 
order to include the motion of the underlying medium we would need to know the full 
solution of the three-dimensional hydrodynamical equations describing the QGP expansion 
and consequently Eo(r, <p, z, t). The appearance of a coordinate dependent quantity in the 
denominator of (1301) would make our expansion of the Euler and continuity equations too 
complicated. A simple way to estimate the effect of the expansion is to represent the cooling 
of the background by the Bjorken formula [2j 



where the proper time is given by r 
(r/t) 2 and thus: 



T(t) 
T(t ) 

t 

7 



70 



1/3 



(55) 



t\/l — v 2 . We have only radial flow v 2 



Vt 2 



(56) 



The initial proper time is taken to be r = 1 fm. With the inclusion of Bjorken cooling the 
term in parenthesis in (1521) will become: 



Tb_ _Tb_ 

To ^ T (r) 



Ty, 



T (r ) 



, r o. 



1/3 



(57) 



Inserting (157]) into wave equation ( |54l) we have: 



dr\ dt s dr 2 



T, 



T (r ) 



To, 



1/3 



\T (r ) \t 

C + 4tA^£i1 
s 3s) dr 2 { 



1/3 



* dii i\ 
£l ~dV + Yt 



c s d 2 i 1 
2 dz 2 



(58) 
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for a viscous fluid. Inserting (1ST)) into ( |52|) we have: 



dr\ dt s dr 2 



1 + 



T, 



B 



T (t q ) \t 



1/3 



£l 97 + 2t 



+ 



2 a? 2 



(59) 



for an ideal fluid. 



V. NUMERICAL RESULTS AND DISCUSSION 



For simplicity, we assume that when they are formed and also throughout the expansion 
the tubes are uniform along the longitudinal direction and therefore: 



c s d 2 ii 







2 dz 2 

Integrating (152]) and ( )54|) with respect to r and setting the integration constant to zero we 
arrive at the cylindrical breaking wave equation for the ideal fluid: 



dii dix c s 
dt +Cs dr + 2 



£l ^ + 2^ = ° 



(60) 



and at the famous Burgers equation 

1+1 



19 



dii dix c s 
dt s dr + 2 



Tb 
To 



42 1| for the viscous fluid: 
l 



(61) 



which in this case is a cylindrical Burgers equation [25|, |26[. Both can be solved numerically 
for a given choice of T and Tb- c 2 s = 1/3 for this equation of state. If both equations 
had only the first two terms, they would describe a traveling wave with velocity c s . The 
third term makes the equations nonlinear in ii. Its effect is to increase the velocity of the 
wave, which is given by the coefficient of the terms proportional to dii/dr. The velocity 
becomes therefore proportional to ii and so the top of the wave travels faster than its 
bottom. Because of this, an initially gaussian pulse turns into a triangular pulse with a 
"vertical wall", as it will be seen in the figures. Finally, the nonlinear term induces rapid 
oscillations around the region close to the wall. This is called dispersion. The term i\/2t 
in both equations comes from the use of cylindrical geometry. It causes the attenuation of 
the initial perturbation at increasing times. Changes in the equation of state imply changes 
in the evolution of the tube. A harder EOS will have a bigger velocity of sound and this 
will make the tube move faster. Moreover the strength of the nonlinear term is directly 
proportional to T B , and consequently (because of ff24|) ) to the bag constant, which, in its 
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turn, contains information about the nonperturbative components of the EOS. Increasing 
the bag constant makes the tube move faster! Inversely, increasing the temperature of the 
background, T , makes the pulse to propagate slower. In spite of the qualitative richness 
of (I6D|) and (I6TT) . for realistic values of B and T , the nonlinear term has a very limited 
range of numerical values. Moreover, as we can observe in ( I6UD and in 0611) . this term is 

on 

never large. Thus we can conclude a posteriori that the linearization, as performed in p, |8J, 
may indeed be a good approximation. In the case of the Burgers equation fl6T|) the second 
order derivative term tames the breaking and dispersion of the wave and at the same time, 
dissipation reduces its amplitude. 

The initial condition is given by a gaussian pulse in e x : 



Ae~ r2/r o (62) 



where the amplitude A and the approximate width r$ are parameters which depend on the 
dynamics of flux tube formation. For simplicity we shall refer to ro as the initial "radius" 
of the tube. If the tubes are perturbations we expect that A < 1. According to current 
estimates 23| the transverse size of the tubes is of the order of 1 fm and thus in our 
calculations 0.1 fm < ro < 0.8 fm. We consider hot QGP at temperatures To = 150 MeV 
and T = 500 MeV treated as an ideal fluid (rj/s = (/s = 0) and as a viscous fluid (rj/s = 



0.08 and (/s = 0) [22|,|23|. 



In the numerical analysis there are many cases to be considered. We present our results in 
eight figures (Figs. 2-9). The first four refer to the static background and the second group 
of four shows solutions for the same parameters for the case of an expanding background. 
All figures have four panels. The two upper panels refer to the low temperature (To = 150 
MeV) and the two lower panels to the high temperature (T = 500 MeV). The two panels 
on the left show results with the ideal fluid and the two panels on the right results with the 
viscous fluid. From the figures we want to see how sensitive the expansion of the tube is to 
changes in: i) the initial amplitude, ii) the temperature, iii) the tube radius, iv) the strength 
of viscosity and v) the expansion of the background fluid. In what follows we discuss the 
role played by each one of these variables mentioning them by order of relevance. 
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1 . Viscosity 



The most striking finding is the strong influence of viscosity. This can be seen in all 
figures and most clearly in the comparison between Fig. [2^) and [2b). Viscosity damps the 
amplitude of the pulse by a factor ten in 1 fm! Increasing the temperature of the medium 
reduces the effect of viscosity as it can be seen from the factor 1/Tq in (!6T|) . However the 
comparison between Fig. [2b) and [2H) shows that this reduction is not very strong. The role 
played by viscosity is also reduced when the initial radius parameter of the tube goes from 
r = 0.1 to 0.8 fm. This is easy to understand looking at f[62l) and then at the right hand 
side of ( 16 ip . A broader initial distribution generates smaller spatial gradients appearing 
in the viscosity term of (|6ip . which becomes smaller. Nevertheless the attenuation of the 
initial pulse remains strong as compared to the ideal fluid case. This situation is illustrated 
in Fig. [31 which is to be compared with Fig. [2j If we increase the initial amplitude from 
A = 0.5 to A = 0.8, the relevance of viscosity remains the same. This can be checked by 
comparing Fig. [2J with Fig. H] and also comparing Fig. [3] with Fig. [51 The introduction of 
the background cooling preserves the difference between ideal and viscous fluids, as can be 
inferred from the comparison between Figs. [2J [31 SI G3 and their analogues with cooling Figs. 
El U\ E and El 

The introduction of viscosity in our calculations is what more strongly changes them. As 
anticipated in the introduction, due to dissipation, viscosity strongly damps and broadens 
the tubes during their expansion and they are more easily mixed with the background fluid, 
loosing their influence on final state particle correlations. This is a robust conclusion of 
our numerical analysis since it remains valid in all situations considered. Moreover viscosity 
prevents the perturbation wave from breaking, as can be seen comparing, for example, Fig. 
[5^) with [5b) or comparing Fig. [7Ji) with [7b). Looking at the time evolution of the peaks of 
the pulses, we can have an idea of the velocity with which they propagate. Comparing the 
left with right side of all figures we can see the velocity of the pulses is only weakly changed 
by viscosity. This velocity is defined by the sound velocity, which in our approach is the 
same both for ideal and viscous fluids. 
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2. Initial radius of the tube 



The solutions of nonlinear differential equations, such as ( 160]) and ( I6T]) . are expected to 
depend on the initial condition. We can check this dependence changing the parameters in 
(I6"2j) and solving again both (16 Op and (1611) . The comparison between Fig. [2] with Fig. [3] and 
also the comparison between Fig. H]with Fig. [5] shows that, after viscosity, changes in the 
initial tube radius are those which most substantially affect the tube evolution. Essentially 
thinner tube are more fragile. They break more easily, developing secondary bumps (called 
"radiation") and/or forming a wall with rapid oscillations at the edge. These instabilities 
may lead to loss of localization and absorption of the tube by the medium. Larger tubes, on 
the other hand, live longer in the plasma. This conclusion may be relevant for the physics 
of particle correlation studied now at RHIC and LHC Moreover, the transverse size of the 
tubes has physica. origins. G.asma flux tubes are typically thinner fl than the tubes 
obtained in event generator based on string models (see [5] and [6J for details). 



3. Initial amplitude of the tube 

Figs. S] and E] are repetitions of [2] and [3J with a larger amplitude. These cases are 
very interesting for us because in perturbations with larger amplitudes the nonlinear effects 
become more important. The reductive perturbation method (RPM) employed here is well 
suited to preserve the nonlinearities of the original equations and transfer them to the 
equations which govern the evolution of perturbations. From the figures we can conclude 
that pulses with larger amplitudes break faster. However the effect is not very pronounced 
because the range of variation considered here is relatively narrow: 0.5 < A < 0.8. 



4- Temperature 

Comparing in all figures the upper panels with the lower panels, we conclude that, in the 
case of ideal fluids, there are only small differences between them. This weak dependence 
on the temperature might have been anticipated from a closer look at the coefficient of the 
nonlinear term in (1601) and in (16 ip . The temperature dependent term in this coefficient 
can vary only between zero and one, changing the overall coefficient at most by a factor 
two (in realistic calculations the range of variation is even narrower because of the limits 
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in the temperature: 150 < To < 500 MeV). This weak sensitivity comes from all the 
developments which led to (160]) and is difficult to say what is more responsible for the 
final result, whether the equation of state or the approximations adopted. In contrast, the 
temperature dependence of the viscosity term in Eq. (I6T1) is slightly stronger. Therefore the 
comparison of upper with lower panels on right side of all figures reveals more pronounced 
differences. In viscous fluids the increase of temperature decreases the amplitude of the 
pulse and makes it live longer. The second derivative term in the Burgers equations does not 
permit the breaking and dispersion of the pulse. However, together with the geometrical term 
ii/2t, it causes the attenuation of the tube. In our calculations the viscosity coefficients (rj 
and £) were kept constant but they may be temperature dependent, enhancing the sensitivity 
of our results to the temperature. 

5. Background expansion 

We have repeated all the calculations replacing (152]) and (1511) by (158]) and (159]) . The 
numerical solution of the latter equations (neglecting the derivatives with respect to z) is 
presented in Figs. 6 to 9. The effect of the Bjorken cooling is to slightly reduce the amplitude 
of the pulse, which can be best seen comparing Fig. 2 with Fig. 6 and Fig. 3 with Fig. 
7. It is a small effect and this is very interesting. The cooling studied here is a crude 
representation of the real three dimensional background fluid expansion. If we had found 
that cooling is important, this would suggest that a realistic treatment of the background 
expansion would change completely the conclusions obtained so far. This seems not to be 
the case. 

A. Final remarks 

If, on one hand, the similarity between the figures is somewhat deceptive (because of the 
weak dependence on the dynamical ingredients), on the other hand they deliver a strong 
message: the tube expands radially with a supersonic velocity and in less than 4 fm/c it 
becomes a "ring", with a hole in the middle. Moreover, by this time the amplitude is al- 
ready reduced by a factor two and the tube (or ring) looses the strength to "push away" the 
surrounding matter. This agrees with the results found in [6] , where the evolution of a tube 
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was studied in a different way. In that work the numerical solution of the hydrodynamical 
equations of the total system (tubular perturbation + background) was obtained, whereas 
here we have isolated the perturbation from the background and written a differential equa- 
tion for it. We have provided an independent check of the results found in [6| with the use 
of a different equation of state. 

An important conclusion of our work is that viscosity strongly affects the propagation of 
perturbations in the quark gluon plasma. This conclusion was obtained with the relativistic 
Navier-Stokes formalism and it would be interesting to check if it remains valid in other 
relativistic theories of viscosity. 



VI. APPENDIX 



The derivation of Eq. (|54|) is quite similar to derivation of Eq. ( )52|) and here we would like 
to add some details. The expansion of Eqs. ( [16]) and ffT8l) in powers of a is straightforward. 
The calculation is faster if we keep terms only up to a 2 and neglect some higher order terms 
even before reaching the identities equivalent to fHUl) . fT4"Tj) and (T4"21) . It is useful to remember 
that: 

7 n £ 1 + (n/2) x v 2 oc 1 + a 2 + ... 

(63) 

The viscous terms in ffTBl) are of the following order in a: 



~ft = 7 3 ^-tt- oc ct 5/2 and V7 = >y 3 vVv oc a 5/2 
dt ' dt ' ' 



rj ?/ (9^(9^7 oc a 4 



r]v l d. 



du^ 



r)v l d IJ 



M dt 

7( — + v ■ V j (ju 



oc a 



(- 

\dt 



oc a 



oc a 



d'y 
~dt 



+ V ■ (7«) 



OC a" 



■qd^^v 1 ) oc a 2 
r/ d^d 1 ^ oc (j 2 



7^ + w- V^) (7«V M 



OC (J 



(64) 
(65) 
(66) 
(67) 
(68) 
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oc a 



~dt 



V ■ (~/v) 



OC cr" 



(69) 



Using ( 1641) to ( 1691) in ( |T6l) and keeping only terms up to 0(a 2 ) we obtain: 







: dp 



+ Tj 



dt 



dt 



+ [C-^v)ld t (V-v) = (70) 



dt 2 / dt 2 

The 7 n (with n = 1,2, 3) factors in (ITU]) are also multiplying u* or its derivative, which 
contributes with at least one power of a. Therefore we consider 7™ = 1 and (ITOl becomes: 



¥ + (*.V)tT 



(e + p) 



+ 



;^+p) 



1 \. 



r/V 2 -y+ (C + -^JV(V-w 



(71) 



The equation above is the simplified version of the relativistic Navier-Stokes equation and 
from ( ITT]) it is easier to derive ( 153]) . 

Analogously we estimate the order (in a) of the viscous terms in ( fl8l) : 

77 / dj \ 2 



1 
T 

T 



T\dtJ 




oc cr 



V7- — (7^) 



a 



(d l U j )djUi oc a 3 



TV3 



$7 



1 2 



+ 7V • v + V7 • V 



oc a 



(72) 

(73) 
(74) 

(75) 
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FIG. 7: The same as of Fig. [3] with background fluid expansion. The curves show numerical 
solutions of ([60]) and (}61~|) with a changing Tq, given by Eq. ([55]) . 
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FIG. 8: The same as of Fig. [5] with background fluid expansion. The curves show numerical 
solutions of ([60]) and ([6T]) with a changing Tq, given by Eq. ([55]) . 
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FIG. 9: The same as of Fig. [5] with background fluid expansion. The curves show numerical 
solutions of ([60]) and ([6T]) with a changing Tq, given by Eq. ([55]) . 
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